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Abstract. In this paper, the concept of generalized spectral function is introduced for finite- 
order tridiagonal symmetric matrices ( Jacobi matrices) with complex entries. The structure 
of the generalized spectral function is described in terms of spectral data consisting of the 
eigenvalues and normalizing numbers of the matrix. The inverse problems from generalized 
spectral function as well as from spectral data are investigated. In this way, a procedure for 
construction of complex tridiagonal matrices having real eigenvalues is obtained. 
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1 Introduction 

Consider the N x N tridiagonal symmetric matrix (Jacobi matrix) with complex entries 
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where for each n, a n and b n are arbitrary complex numbers such that a n is different from zero: 



Q"n ) bn £ 



a n / 0. 



;i.2) 



In the real case 

a n ,b n GR, a n ^ 0, (1.3) 

the matrix J is Hermitian (self-adjoint) and in this case many versions of the inverse spectral 
problem for J have been investigated in the literature, see [U [2j [3] and references given therein. 

In the complex case fjl .2[) . the matrix J is in general non-Hermitian (non-selfadjoint) and our 
aim in this paper is to introduce appropriate spectral data for such a matrix and then consider 
the inverse spectral problem consisting in determining of the matrix from its spectral data. 

*This paper is a contribution to the Proceedings of the Vllth Workshop "Quantum Physics with Non- 
Hermitian Operators" (June 29 - July 11, 2008, Benasque, Spain). The full collection is available at 
|http:// ww.emis.de/journals/SIGMA/PHHQP2008.htmf1 
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As is known [U [6j [9] , for non-selfadjoint differential and difference operators a natural 
spectral characteristic is the so-called generalized spectral function which is a linear continuous 
functional on an appropriate linear topological space. In general very little is known about the 
structure of generalized spectral functions. 

Given the matrix J of the form (jl.ip with the entries satisfying (|1.2p . consider the eigenvalue 

=0 ' 

difference equation 



problem Jy = Xy for a column vector y = {y n }n=o^ that is equivalent to the second order linear 



a n _iy n _i + b n y n + a n y n+1 = Xy n , n G {0, 1, . . . , N - 1}, a_i = a N -i = 1, (1.4) 
for {y n }n=-ii with the boundary conditions 

y-i = vn = o. (1.5) 

The problem (II. 4h . (jl,5p is a discrete analogue of the continuous eigenvalue value problem 

+ q(x)y(x) = Xy(x), x e [a,b], (1.6) 



p{x)-^-y{x) 
ax 



d_ 

dx 

y(a)=y(b)=0, (1.7) 



where [a, 6] is a finite interval. 
To the continuous problem 



d_ 

dx 
2/(0) = 0, 



P(x)^y(x) 



+ q{x)y{x) = Xy(x), x £ [0, oo), 



on the semi-infinite interval [0, oo) there corresponds a Jacobi matrix of the type (jl.ip but J 
being infinite both downwards and to the right. To the equation in (11.61) considered on the whole 
real axis (— oo, oo) there corresponds a Jacobi matrix which is infinite in the all four directions: 
upwards, downwards, to the left, and to the right. 

The case of infinite Jacobi matrices was considered earlier in the papers [6j [7J [9] in which 
the generalized spectral function was introduced and the inverse problem from the generalized 
spectral function was studied. However, in the case of infinite Jacobi matrices the structure of 
the generalized spectral function does not allow any explicit description because of complexity 
of the structure. Our main achievement in the present paper is that we describe explicitly the 
structure of the generalized spectral function for the finite order Jacobi matrices (jl.ip . (|1.2p . 

The paper is organized as follows. In Section[2l the generalized spectral function is introduced 
for Jacobi matrices of the form (jl.ip with the entries satisfying (jl.2p . In Section [3l the inverse 
problem from the generalized spectral function is investigated. It turns out that the matrix (jl.ip 
is not uniquely restored from the generalized spectral function. There are precisely 2 N_1 dis- 
tinct Jacobi matrices possessing the same generalized spectral function. The inverse problem 
is solved uniquely from the data consisting of the generalized spectral function and a sequence 
{o"i,o"2, • • • , cr^v-i} of signs + and — . Section [4] is devoted to some examples. In Section [5j we 
describe the structure of the generalized spectral function and in this way we define the concept 
of spectral data for matrices (jl.ip . In Section [U the inverse problem from the spectral data is 
considered. In Section [7J we characterize generalized spectral functions of real Jacobi matrices 
among the generalized spectral functions of complex Jacobi matrices. In Section [H we describe 
the structure of generalized spectral functions and spectral data of real Jacobi matrices. Finally, 
in Section [91 we consider inverse problem for real Jacobi matrices from the spectral data. 

Note that considerations of complex (non-Hermitian) Hamiltonians in quantum mechanics 
and complex discrete models have recently received a lot of attention [10j EU EH Q3] . For some 
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recent papers dealing with the spectral theory of difference (and differential) operators with 
complex coefficients see [HI E51 ESI E3 CES]- Fur further reading on the spectral theory of the 
Jacobi difference equation (three term recurrence relation) the books [191 [201 [21] [22l [23] are 
excellent sources. 



2 Generalized spectral function 

Given a matrix J of the form (jl.ip with the entries satisfying (jl.2p . Consider the eigenvalue 
problem Jy = Xy for a column vector y = {y n } n =o, that is equivalent to the second order linear 
difference equation 

a n -iy n -i + b n y n + a n y n+ i = \y n , n e {0, 1, . . . , N - 1}, o_i = a N _i = 1, (2.1) 
for {y n }n=-ii with the boundary conditions 

y_i = yv = 0. (2.2) 
Denote by {i-'„(A)}^ r = „ 1 the solution of equation (|2.ip satisfying the initial conditions 

y_i = 0, y = 1. (2.3) 

Using (|2.3p . we can find from equation (|2.1I) recurrently the quantities -P n (A) for n = 1, 2, . . . , N; 
P n (X) is a polynomial in A of degree n. 

Thus {P n {X)}n=o is the unique solution of the recursion relations 



b P (X)+a P 1 (X) = XP (X), 

a„_iP n _i(A) + b n P n (X) + a n P n+ i(A) = AP n (A), 

n G {1,2,..., JV- 1}, a N -i = l, 



(2.4) 



subject to the initial condition 

Po(X) = 1. 
Lemma 1. T/ie equality 

det (J - A/) = (-lj^aoai • • • a Ar _ 2 P 7V (A) 
holds so that the eigenvalues of the matrix J coincide with the zeros of the polynomial Pn(X) 



(2.5) 
(2.6) 
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and set A n (A) = det(J n — XI). Note that by / we denote an identity matrix of needed size. By 
expanding the determinant det(J n+ i — XI) by the elements of the last row, it is easy to show 
that 



A n+ i(A) = (b n - A)A„(A) - al^A^tiX), 



n 



1,2,... 



A (A) = 1. 
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Dividing this equation by the product ao • • • £t re -i> we find that the sequence 

d-i = 0, d = 1, d n = (-l)"(a • • • a„_i)" 1 A n (A), n = 1, 2, . . . , 

satisfies (|2.ip . f|2.3jl . Then d„ = -P„(A), re = 0, 1, ... , and hence we have (|2.6|) because Jjy = J 
and Q TV— l = 1- B 

For any nonnegative integer m denote by C m [A] the ring of all polynomials in A of degree 
< m with complex coefficients. A mapping : C m [A] — > C is called a linear functional if for 
any G(X), H(X) G C m [A] and a G C, we have 

(n, G(A) + #(A)) = (fi, G(A)> + (J2, #(A)) and (fi, aG(X)) = a (fi, G(A)} , 

where (0,G(A)) denotes the value of Q on the element (polynomial) G(A). 

Theorem 1. There exists a unique linear functional £2 : C2n[X] — > C suc/i i/iaf i/ie relations 

(n,P m (X)P n (X)) = 5 mn , m,n€{0,l,...,JV-l}, (2.7) 
(n,P m (A)Pjv(A)) = 0, m€{0,l,...,JV}, (2.8) 

ZioW, where 5 mn is the Kronecker delta. 

Proof. First we prove the uniqueness of f2. Assume that there exists a linear functional Q 
possessing the properties (|2.7[) and f|2.8j) . The 2 A + 1 polynomials 



P n (X) (n = 0,l,..-,A-l), P m (X)P N (X) (m = 0,l,...,A) (2.9) 

form a basis for the linear space C2/v[A] because they are linearly independent (their degrees 
are distinct) and their number 2 A + 1 = dimC2/v[A]. On the other hand, by (|2.7p and (|2.8p the 
functional takes on polynomials (|2.9p completely definite values: 



<n,P„(A)) = 5on, nG{0,l,...,A-l}, (2.10) 
(n,P m (X)P N (X)) =0, mG {0,1,..., A}. (2.11) 

Therefore O is determined on C2/v[A] uniquely by the linearity. 

To prove existence of f2 we define it on the basis polynomials (|2.9p by (|2.10p . (|2.1ip and then 
we extend to over the whole space C2tv[A] by linearity. Let us show that the functional Q 
defined in this way satisfies (12. 7b . (12. 8p . Denote 



(n,P m (X)P n (X))=A mn , m,ne {0,1,..., A}. (2.12) 

Obviously, A mn = A nm for m, n G {0, 1, . . . , A}. From (|2.10p and (|2.1ip we have 

A m0 = A 0m = 5 m0 , me{0,l,...,4 (2.13) 

A mN = A Nm = 0, mG {0,1,..., A}. (2.14) 

Since {Pti{X)}q is the solution of equations (|2.4p . we find from the first equation, recalling (|2.5p . 

X = b + a Pi(X). 

Substituting this in the remaining equations of (|2.4|) . we obtain 

a n _iP n _i(A) + b n P n {X) + a n P n+1 (X) = b P n (X) + a Pi(X)P n (X), 
n G {1,2,..., A- 1}, ajv -i = l. 
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Applying the functional to both sides of the last equation, and recalling (|2.13p and (|2.14p . we 

get 

A nl = A ln = 5 nU ne{0,l,...,JV}. (2.15) 
Further, since 

bmPm (A) + (X) = XP m (X), me {1,2,...,JV-1}, 

a n -iP n ~i(X) + b n P n (X) + a n P n+ i(X) = AP n (A), n € {1, 2, . . . , JV - 1}, 

we obtain, multiplying the first of these identities by P n (X), and multiplying the second by 
P m (X), then subtracting the second result from the first: 

a m -iPm-i{X) P n (X) + b m P m (X)P n (X) + a m P m+ i(X)P n (X) 

= a n _iP n _!(A)P m (A) + b n P n {X)P m {X) + a n P n+ i(A)P m ,(A), m,ne{l, 2, ...,N-1}. 

Applying the functional f2 to both sides of the last equation, and recalling (|2.13p . (|2.14p . 
and (|2.15p . we obtain for A mn the boundary value problem 

"m-l-^m-l,n b m Amn ~\~ Q"mAm+l,n — ^n— \A n — ~\~ b n A nm -\- Q, n ^4 n ^i j?Tl , (2.16) 

m,7iG{l,2,...,JV-l}, 
A n Q = A 0n = 5 n Q, Ani = A\ n = 8 n \, A^n = A n N = 0, (2-17) 
nG{0,l,...,iV}. 

Using (|2.17|) . we can find from (|2.16|) recurrently all the A mn and the unique solution of problem 
(pnE]) . (pTTP is Arm = S mn for m, n e {0, 1, . . . , N- 1} and A mN = for m € {0, 1, ... , N}. U 

Definition 1. The linear functional 17 defined in Theorem [T] we call the generalized spectral 
function of the matrix J given in (jl.ip . 

3 Inverse problem from the generalized spectral function 

The inverse problem is stated as follows: 

1. To see if it is possible to reconstruct the matrix J, given its generalized spectral function 0. 
If it is possible, to describe the reconstruction procedure. 

2. To find the necessary and sufficient conditions for a given linear functional 0, on C2iv[A], 
to be the generalized spectral function for some matrix J of the form (|l.ip with entries 
belonging to the class (jl.2p . 



Since P n (X) is a polynomial of degree n, we can write the representation 



n-l 



Pn(X) = a n [X n + Y,XnkX k ) , TIG {0,1,..., JV}. (3.1) 



fc=0 



Substituting (|3,ip in (|2.4p . we find that the coefficients a n , b n of system (j2.4|) and the quanti- 
ties a n , Xnk of decomposition (|3.1|) . are interconnected by the equations 

QL 

(0<ra<iV-2), ao = l, a N = a N - X , (3.2) 



«n+l 

frn = Xn,n-1 ~ Xn+l,n (0 < n < JV - 1), X0,-1 = 0. (3.3) 
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It is easily seen that relations (|2.7p , (|2.8p are equivalent to the collection of the relations 

(n,X m P n (X)} = m = 0,l,...,n, tiG{0,1,...,JV-1}, (3.4) 

(fl,X m P N (X)) = 0, m = 0,1,... ,iV. (3.5) 
In fact, using (|3.ip . we have 

m— 1 

(fi, P m (A)P n (A)) = a m (0, A m P n (A)> + a m ]T Xmj («, A J P„(A)> . (3.6) 

3=0 

Next, since we have the expansion 

j 



i=0 



\i = Y / c<f>P i (X), jG{0,l,...,iV}, 



it follows from (|3.6p that (|3.4p . (|3.5p hold if we have (|2.7p . (j2.8|) . The converse is also true: 
if ([3~4"|) . (1331) hold, then ([277]), can be obtained from ([32]), in conjunction with (pTTjh 
Let us set 

S/ = (0,A'), Ze{0,l,...,2iV}, (3.7) 

that are the "power moments" of the functional 0. 

Replacing P n (X) and Pjv(A) in (|3.4|) and (|3.5|) by their expansions in (|3.ip . we obtain 

n-l 

s n+m + ^ XnkSk+m = 0, m = 0,1,..., n-l, n G {1,2, . . . , N}, (3.8) 

fc=0 
7V-1 

S2N + ^ XNkSk+N = 0, (3.9) 
fc=0 

n— 1 ^ 

S2n + y^XnkSk+n = —o, n G {0, 1, . . . , JV - 1}. (3.10) 

k=0 n 

Notice that (|3.8p is the fundamental equation of the inverse problem, in the sense that it 
enables the problem to be formally solved. For, if we are given the linear functional on C2tv[A], 
we can find the quantities S[ from (|3.7p and then we consider the inhomogeneous system of linear 
algebraic equations (|3.8p with unknowns Xn0> Xni, ■ ■ ■ , Xn,n-i> f° r every fixed n€ {1,2,..., N}. 

n— 1 

If this system is uniquely solvable, and S2 n + ^2 Xnk s k+n 7^ for n G {1, 2, . . . , N — 1}, then the 

fc=o 

entries a n , of the required matrix J can be found from (|3.2p and (|3.3p . respectively, a n being 
found from (|3,10p . The next theorem gives the conditions under which the indicated procedure 
of solving the inverse problem is rigorously justified. 

Theorem 2. In order for a given linear functional VL, defined on C2iv[A], to be the generali- 
zed spectral function for some Jacobi matrix J of the form (II. ip with entries belonging to the 
class (jl.2p . it is necessary and sufficient that the following conditions be satisfied: 



(i) (fi, 1) = 1 (normalization condition); 
(ii) if, for some polynomial G(X), degG(A) < N — 1, 

(n,G(X)H(X)) = 



for all polynomials H(X), degH(X) = degG(A), then G(X) = 0; 
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(Hi) there exists a polynomial T(X) of degree N such that 
(n,G(X)T(X)) = 
for all polynomials G(X) with degG(X) < N. 

Proof. Necessity. We obtain (i) from (|2.7p with n = m = 0, recalling (|2.5p . To prove (ii), we 
write the expansion 

m 

G(X) =^4 m) P i (A), m = degG(A), 
i=o 

and take as H(X) the polynomial 

m 

H(x) = ^c { r ) m), 

3=0 

where the bar over a complex number denotes the complex conjugation. Then we find from 
(n,G(X)H{X)) = using (JUT]) that 



3=0 



I M 2 



hence c[ m ^ = 0, j = 0,1,..., m, i.e., G(A) = 0. The statement (Hi) of the theorem follows 
from ([22|) if we take T(A) = Pat (A). 

Sufficiency. The proof will be given in several stages. 

(a) Given the linear functional 0, defined on C2iv[A] and satisfying the conditions of the 
theorem. Consider equation (|3.8p with the unknowns Xnk, k = 0, 1, ... ,n — 1, in which s; are 
found with the aid of the functional f2 from expression ()3.7I) . Let us show that this equation 
has a unique solution for every fixed n £ {1,2, ... , N}. For this, it is sufficient to show that the 
corresponding homogeneous equation 



n-l 



^2 9kSk+m = 0, m = 0,1, ... ,n- 1, (3.11) 

fc=0 

has only the zero solution for every n. Assume the contrary. For some n £ {1,2, . . . ,N} let 
equation (|3.1ip have the nonzero solution (gk)^ 1 ■ Further let (/i m )o _1 be an arbitrary vector. 
We multiply both sides of (I3.11|) by h m and sum over m between and n — 1; we get 

n— 1 n— 1 

^ ^ hm9kSk+m = 0. 
m=0 fc=0 

Substituting expression (13.71) for Sfc +m in this equation and denoting 

n— 1 n— 1 

G(A) = ^ 5fc A fc , (A) = J> OT A™, 

fc=0 m=0 

we obtain 



(n,G(X)H(X)) = 0. 



(3.12) 
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Since (/im)^ 1 is an arbitrary vector, we find from (|3.12p . in the light of condition (ii) of the 
theorem, that G(X) = 0, and hence go = g\ = ■ ■ ■ = g n -i = 0, in spite of our assumption. Thus, 
for any n £ {1,2, ... , N}, equation (|3.8p has a unique solution, 
(b) Let us show that 



n-1 



S2n + ^2xnkSk+n^0i n € {1,2, . . . ,N - 1}, (3.13) 



k=0 

where (Xnk^Zo * s ^ e solution of the fundamental equation f|3.8j) . (For n = 0, the left-hand side 
of (|3.13p is sq = (^j 1) = !•) Assume the contrary, i.e., for some n E {1, 2, . . . , N — 1} 

n-1 

s 2n + XnkSk+n = 0. 
fc=0 

Joining this equation to the fundamental equation (|3.8|) . we obtain 

n-1 

s n+m + ^ XnkSk+m = 0, m = 0, l,...,n. (3.14) 

fc=o 

Let (/i m )o be an arbitrary vector. Multiplying both sides of (I3.14p by h m and summing over m 
from to n, we obtain 

n n n—1 

m=0 m=0 fc=0 

Replacing s/ in this by its expression (|3.7p . we obtain 
(n,[A n + X (A)]H(A)) = 0, 

where 

n—1 n 

x(A) = ]>> nfe A fc , iT(A) = £ ^A m . 

k=0 m=0 

Since (/i m )o is an arbitrary vector, we obtain from the last equation in the light of condition (ii) 
of the theorem: 

A" + x(A)^0, 

which is impossible. Our assumption is therefore false. 

(c) Given the solution (Xnfc)fc=o °f the fundamental equation (|3.8p . we find a n from (|3.10p 
for n G {0, 1, . . . , N — 1} with cto = 1 and set a^v = ckjv-i- Then we find the polynomials -P n (A) 
from j33|). Let us show that the relations flUT]), <^8^ hold. It is enough to show that ([33]), §2M 
hold because (^4l) . (|3T5|) together are equivalent to relations (|2~7l) . (p^l) . From (^8|) and (i3~T0l) 
we have (|3.4p and ()3.5|) . the latter except for m = N. So it remains to show that 

(fi, X N P N (X)) = 0. 

For this purpose we use the condition (m) of the theorem. By this condition we have 

N 



k=0 
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and 

N 

= (n,p N (x)T(\)) = J2t k (n,\ k p N (X)) =t N (n,x N p N (\)), 

k=0 

where we have used (|3.5p except for m = N. Hence (Q, X N Pat (A)) = 0. 

(d) Let us show that the polynomials P n (A), n = 0,1,..., N, constructed in accordance 
with (|3.ip with the aid of the numbers Xnk and a n obtained above, satisfy the equations 

6 Po(A) + a Pi(A) = AP (A), 

a n _iP n _i(A) + b n P n (X) + a n P n+l {X) = XP n (X), (3.15) 
n G {1,2,..., iV- 1}, ajv-i = l, 

where the coefficients a n , b n are given by the expressions 

OL 

a n = — — (0<n<N — 2), ao = 1, ajv = a/v-i, (3.16) 

K = Xn,n-1 ~ Xn+l,n (0 < 71 < N - 1), X0 -1 = 0. (3.17) 

We first verify the first equation of (|3.15p . From (|3.ip we have 

P (A) = 1, Pi (A) = ai(A + X io). 
Hence the first equation of (|3.15p has the form 

bo + a ai(X + xio) = A; 
and this is true, since, by (|3.16p and (|3.17p . 

a ai = 1, b = -xio- 

Let us prove the remaining equations of (|3.15p . Since AP n (A) is a polynomial of degree n + 1, 
while Pfe(A), k = 0, 1, . . . , n + 1, are linearly independent, we have 

n+l 

AP n (A) = ^4 n) P fc (A) ne{l,2,...,N-l}, c^ _1) = 1, (3.18) 

fc=0 

where c£ , k = 0, 1, . . . ,n + 1, are constants. By (|2.7p . (|2.8p which we proved in (c), we have 
from (1338]) : 

4 n) = (n,AP n (A)Pfc(A)), fc = 0,l,...,n + l (n G {1,2, ... ,7V — 2}). (3.19) 

The polynomials APfc(A), k = 0, 1, . . . , n — 2, have degrees < n— 1, and hence we find from (13.191) 
in the light of (JUT]), ([HI]) that 

4 n) = 0, fc = 0,l,...,n-2 (ne{l,2,...,iV-l}). 
Consequently, expansion (|3.18p takes the form 

c^Pn-^X) +cWp n (A) + C i n iP n+1 (A) = AP n (A), 7i £ {1, 2, . . . , N — 1}. (3.20) 
It follows from (|3. 19j) that c^-i = &n ■ Hence, denoting 

ci n i=Sn, cW=6 n , (3.21) 
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we have from (|3.20[) : 

S„_iP„_i(A) + b n P n {\) + a n P n+1 {\) = XP n (X), 



n £ {1,2,...,N - 1}. 



(3.22) 



Replacing P n (X) in (|3.22|) by its expression (|3.ip and equating coefficients of like A n , while 
recalling (|3.16p . (|3.17p . we obtain 



a. 



a,, 



bn = Xn,n— 1 ~~ Xn+1 



a n (0<n<iV-2), 

b n (0<n<N -1). 



Theorem [2] is completely proved. ■ 

Remark 1. It follows from the above solution of the inverse problem that the matrix (jl.ip is 
not uniquely restored from the generalized spectral function. This is linked with the fact that 
the a n are determined from (13.10p uniquely up to a sign. To ensure that the inverse problem 
is uniquely solvable, we have to specify additionally a sequence of signs + and — . Namely, let 
{a\, o~ 2 , . . . , cttv-i} be a given finite sequence, where for each n & {1, 2, . . . , N — 1} the a n is + 
or — . We have 2 N ~ l such different sequences. Now to determine a n uniquely from (|3.10p for 
n € {1, 2, . . . , N — 1} (remember that we always take «o = 1) we can choose the sign o~ n when 
extracting the square root. In this way we get precisely 2 N ~ l distinct Jacobi matrices possessing 
the same generalized spectral function. For example, the two matrices 



" 1 1 " 




1 -1 " 


1 1 




-1 1 



as well as the four matrices 



1 


1 







1 




1 






-1 















-1 

1 1 
1 1 



1 1 
1 1 

-1 








1 -1 


-1 




-l l -l 


1 




-1 1 



have the same generalized spectral function. The inverse problem is solved uniquely from the 
data consisting of Q and a sequence {a\, oi, ■ ■ ■ , otv-i} of signs + and — . 



Using the numbers 

si = (n,\ l ), z = o,i,...,2iv, 

let us introduce the determinants 



(3.23) 



D r , 



so 
si 



•51 
•52 



Sn+1 



n 



0,1,... ,N. 



(3.24) 



It turns out that Theorem [2] is equivalent to the following theorem. 

Theorem 3. In order for a given linear functional £1, defined on C2at[A], to be the generali- 
zed spectral function for some Jacobi matrix J of the form (jl.lj) with entries belonging to the 
class (jl.2p . it is necessary and sufficient that 



D = l, D n ^0 (n = l,2,...,iV-l), 
where D n is defined by (|3.24p and ()3.23j) . 



and 



D 



N 



0. 



(3.25) 
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Proof. Necessity. The condition Dq = 1 follows from 1 = (0, 1) = sq 
for a polynomial 

G(\)=Y,9k\ k 

fc=0 

with n < N — 1 we have 

(n,G(X)H(X)) = 

for all polynomials 
ii 

H(X) = £ h m X m , 

m=0 

then G(X) = 0, that is, go = g\ = • • • = g n = 0. 

If we substitute (13361) and (ET281) in (I3T271) . then we get 



Z>o- By Theorem [2J if 
(3.26) 

(3.27) 
(3.28) 



I 9kSk+m 
\fc=0 > 



0. 



m=0 



Since ho,h\, . . . ,h n are arbitrary, the last equation gives 

n 

^5fcSfc +m = 0, m = 0,l,...,ra. 



(3.29) 



fc=0 



This is a linear homogeneous system of algebraic equations with respect to go, g±, . . . ,g n and the 
determinant of this system coincides with the determinant D n . Since this system has only the 
trivial solution go = g\ = ■ ■ ■ = g n = 0, we have that D n ^ 0, where n < N — 1. 
To prove that Dn = 0, we write equation (|3.8|) for n = N to get 



N-l 



s N+m + XNkSk+m = 0, m = 0,1, ... ,JV- 1. 

, . . . , XN,N-i- Next, these equalities together 



k=0 



This equation has the unique solution xno,Xni 
with (13.91) can be written in the form 



SJV 




so 




Sl 


S7V+1 




Sl 




S2 




+ X^O 




+ XN1 




S2JV-1 




SN-l 




SN 


S2N 




S N 




SN+1 



+ 



Xn,n-i 



sn-i 

S2N-2 
S2N-1 

This means that the last column in the determinant Dn is a linear combination of the remaining 
columns. Therefore Dn = 0. 

Sufficiency. Given the linear functional Q : C2tv[A] — > C satisfying the conditions (|3.25 j> . it 
is enough to show that then the conditions of Theorem [2] are satisfied. We have (fi, 1) = so = 
Do = 1. Next, let (|3.27p hold for a polynomial G(X) of the form (|3.26p and all polynomials H(X) 
of the form (|3.28|) . Then (|3.29|) holds. Since the determinant of this system is D n and D n ^ 
for n < N — 1, we get that go = gi = • • • = g n = 0, that is, G(X) = 0. Finally, we have to show 
that there is a polynomial T(A) of degree N such that 



(n.G(A)r(A)) = o 



(3.30) 
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for all polynomials G(A) with degG(A) < N. For this purpose we consider the homogeneous 
system 



N 



^2t k s k+m = 0, m = 0,1, . . . ,N, (3.31) 



k=0 



with the unknowns to, t±, . . . ,tjv- The determinant of this system is D^. Since by condition 
Dn = 0, this system has a nontrivial solution to, t±, . . . , ijv- We have tw 7^ 0. Indeed, if tjy = 0, 
then we get from (|3.31|) 



N-l 



Y^tkS k+m = 0, m = 0,l,...,N-l. (3.32) 



fc=0 



The determinant of this system is -Djv_i and by condition Z?at-i ^ 0. Then to = t± = ■ ■ ■ = 
t = and we get that the solution to,t±, . . . ,tw of system f|3.31 j) is trivial, which is a contra- 
diction. Taking the nontrivial solution to, t±, . . . , ijv of system (13.3ip we construct the polynomial 



N 



k=0 

of degree N. Then substituting Sk+ m = {£l,\ k+m ) in (|3.3ip gives 

(n,\ m T(\)) = 0, m = 0,l,...,N. 
Hence (|3~30jl holds for all polynomials G(A) with deg G(A) < N. ■ 

Note that the determinant of system (13. 8p coincides with D n -\. Denote by (k = 

0,1,..., m) the determinant that is obtained from the determinant D m by replacing in D m 
the (k + l)th column by the column with the components s m +i, s m+ 2, ■ ■ ■ , S2m+i- Then, solving 
system (|3,8p by making use of the Cramer's rule, we find 

D n-1 

Xnk = — ft~> = 0, 1, . . . ,n - 1. (3.33) 

Next, substituting the expression (|3.33p of Xnk into the left-hand side of (|3.10p . we get 

a- 2 = D n D- l _ v (3.34) 

Now if we set = A m , then we get from (I3T61) . (I3T7D . by virtue of (l3T33|) . (13341 . 

a n = ±( J D n _i J D n+ i) 1/2 D- 1 , nG {0, 1, . . . , iV - 2}, D_i = 1, (3.35) 
6 n = AnD- 1 - A^Zr^, ne {0, 1, . . . ,iV - 1}, A_i = 0, A = si- (3.36) 

Thus, if the conditions of Theorem [3] or, equivalently, the conditions of Theorem [2] are sa- 
tisfied, then the entries a n , b n of the matrix J for which Q is the generalized spectral function, 
are recovered by the formulas (|3.35p . (|3.36p . where D n is defined by (|3.24p and (|3.23p . and A n 
is the determinant obtained from the determinant D n by replacing in D n the last column by 
the column with the components s n +i, s n +2> ■ ■ ■ , S2n+i- 
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4 Examples 

In this section we consider some simple examples to illustrate the solving of the inverse problem 
given above in Section [3j 

Example 1. The functional 
(n,G(A)> = f G(X)dX 



o 



satisfies the conditions (i) and (ii) of Theorem [2j but it does not satisfy the condition (Hi) of 
this theorem. 

In fact, obviously, (CI, 1) = 1. Next, let for a polynomial 
N-l 

G(X) = 9kX k (4.1) 

fc=0 

we have 



(Cl, G(X)H(X)) = [ G(X)H(X)dX = 
Jo 



for all polynomials 

N-l 



H(X) = h k X k , degH(X) = degG(A). (4.2) 

Taking, in particular, 
N-l 

H(X) = Y9k* k , (4-3) 

fc=0 

where the bar over a complex number denotes the complex conjugation, we get 
\G{X)\ 2 dX = 

o 

and hence G(X) = 0. 

The same reasoning shows that there is no nonidentically zero polynomial T(X) such that 
(n,G(X)T(X)) = for all polynomials G(X) with degG(A) < degH(X). 

Example 2. The functional 

JV 



(Q,G(X)} = Yc k G(X k ) 



k=l 

where Ai, . . . , Av are distinct real numbers, c\, . . . ,cjv are complex numbers such that 

JV 



Ck = 1 and Re Cfc > (k = 1, . . . , N), 



k=l 

satisfies the conditions of Theorem [2j 
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In fact, obviously, (Q, 1) = 1. Next, assume that for a polynomial G(A) of the form (|4.ip we 
have 



JV 



<fi, G(X)H(X)) = c k G(X k )H(X k ) = 



k=l 



for all polynomials H(X) of the form (|4,2p . If we take, in particular, H(X) of the form (|4.3p . 
then we get 



v 



^c fe |G(A fe )| 2 = 0. 



k=l 

Hence, taking the real part and using the condition Re c& > (k = 1, . . . , N), we get G(Afe) = 0, 
k = 1, . . . , N. Therefore G(X) = because Ai, . . . , Xn are distinct and G(X) is a polynomial with 
degG(A) < JV - 1. 

Further, for the polynomial 

T(A) = (A-Ai)---(A-Ajv) 

we have (Q, G(X)T(X)) = for all polynomials G(X) so that the condition (Hi) of Theorem [2] is 
also satisfied. Thus the functional satisfies all the conditions of Theorem [2J 
Consider the case TV = 2 and take the functional f2 defined by the formula 

(fi, 67(A)) = cG(0) + (1 - c)G(l), 

where c is any complex number such that c 7^ and c / 1. Let us solve the inverse problem for 
this functional by using formulas (|3.35p and (13.36p . We have 

1,2,..., 



so = 


(n,i) 




1 




Sl 


= (n,x l ) 




1 - 


- c 


for all Z 






1, 




D = 


so 


= 1 














Di = 




so 

Si 


si 

S2 




1 

1 - c 


1 - 
1 - 


c 
c 




c(l 


-c) 








so 


Si 




S2 






1 


1 - 


- c 


1 - 


- c 




D 2 = 




Si 


S2 


S3 




1 


— c 


1 - 


- c 


1 - 


- c 


= 0, 






S2 


«3 


S4 




1 


— c 


1 - 


- c 


1 - 


- c 




A_i 




0, 




H 


^0 = 


Sl 


= 1 


- c, 












Ai = 






so 

Sl 


S2 
S3 




1 

1 - 


c 


h- 1 1— 1 


- c 

- c 


= c(l-c). 



Therefore the functional satisfies all the conditions of Theorem[3l According to formulas (|3.35l) 
and (I3.36|) . we find 



a = ±{D_ 1 D 1 ) 1 / 2 Dq 1 
bo 



Ao^o" 1 



= 1 



- c, 

1 - (1 - c) = c. 

Therefore there are two matrices J± for which 0, is the spectral function: 



Ai^r 1 



AoDq" 1 



" 60 


a 




a 







1 



±v/c(l-c) 



cj c 

The characteristic polynomials of the matrices J± have the form 
det(J± - XI) = A(A - 1). 
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Example 3. Let N = 2. Consider the functional VI defined by the formula 
(n,G(X)) = G(Xo) + cG'(X ), 

where Ao and c are arbitrary complex numbers such that c / 0. This functional satisfies all the 
conditions of Theorem[2l As the polynomial T(X) presented in the condition (Hi) of Theorem[2l 
we can take T(A) = (A — Aq) 2 - 



We have 



so 


= (0,1) 




1 




si -- 


£>_ 


l - 


1, 




D = 


so = 


£>i 




so 

Sl 


Sl 
S2 




Ao 






so 


Sl 




S2 




D 2 




Sl 
S2 


S2 
«3 


S3 
S4 




A_ 


l = 


o, 




A = 


Sl = 


Ai 


= D? 




so 

Sl 


S2 
S3 



(n,x l ) 



A + cIAq" 1 



for 1 = 1,2,..., 



An + c 



1 



A + c A 2 , + 2cA 
A + c Ag + 2cA Ajj + 3cXl 
Ag + 2cA X 3 + 3cXl Ag + 4cAg 

A + c, 

1 



0. 



A^ + 2cA 
A + c Ajj + 3cAg 



2c 2 A . 



Therefore the functional f2 satisfies all the conditions of Theorem[3l According to formulas (|3.35D 
and (|3.36p . we find 



±(XU£>i) 1/2 £>o 1 = ±\fD\ = iV^c 2 = ±ic, 



Oo - 

3 = ^O^o" 1 

h = Ax^f 1 - AoDq" 1 



- A_iDZ 1 1 = A + c, 
-2c 2 An 



(A + c) = A 



c. 



Therefore the two matrices J± for which O is the spectral function have the form 
J± 



' b 


ao 




Ao + c ±ic 








±ic Aq — c 



The characteristic polynomials of the matrices J± have the form 

det(J± - XI) = (A - A ) 2 . 

Note that if N = 3, then the functional 

<n, G(X)) = G(X ) + ciG'(Ao) + c 2 G"(Ao), 

where Ao, c\, c 2 are complex numbers, satisfies the conditions of Theorem [2] (or Theorem [3]) if 
and only if 

c 2 + 0, 2c 2 - c\ £ 0. 



5 Structure of the generalized spectral function 
and spectral data 

Let J be a Jacobi matrix of the form (jl.ip with the entries satisfying (|1.2p . Next, let £1 be the 
generalized spectral function of J, defined above in Section [2j The following theorem describes 
the structure of O. 
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Theorem 4. Let Xi,...,X p be all the distinct eigenvalues of the matrix J and mi, . . . , m p be 
their multiplicities, respectively, as roots of the characteristic polynomial (|2.6|) . There exist 
complex numbers (3k j (j = 1, • • • , m*;, fc = 1, . . . ,p) uniquely determined by the matrix J such 
that for any polynomial G(X) £ C2tv[A] the formula 

P "i k 

(WW) = EEt731v g0 '" 1) ( a *)' (5 - 1} 
k=ij=i u h 

holds, where G^ n \X) denotes the nth order derivative of G(X) with respect to X. 

Proof. Let J be a matrix of the form (jl.ip . (jl.2p . Consider the second order linear difference 
equation 

a n _iy n _i + b n y n + a n y n+1 = Xy n , n € {0, 1, . . . , JV - 1}, a_i = ojv-i = 1, (5.2) 

where {y n }^L_i is a desired solution. Denote by {P n (A)}^ r __ 1 and {Q n (A)}^L_ 1 the solutions 
of equation (15, 2h satisfying the initial conditions 



P_i(A) = 0, Po(A) = l; (5.3) 
Q_i(A) = -l, Q (A) = 0. (5.4) 

For each n > 0, P n (X) is a polynomial of degree n and is called a polynomial of first kind (note 
that -Pn(A) is the same polynomial as in Section [2]) and Q n (X) is a polynomial of degree n — 1 
and is known as a polynomial of second kind. 
Let us set 

M(A) = -f^l. (5.5) 

Then it is straightforward to verify that the entries R nm {X) of the matrix R{X) = ( J — A/) -1 
(resolvent of J) are of the form 

R (\) = l P «( A )^-( A ) + M(X)P m (X)], 0<n<m<N-l, 
nm[ ' \ P m (X)[Qn(X) + M(X)P n (X)}, 0<m<n<N-l. l ' j 

Let / be an arbitrary element (column vector) of C N , with the components /o, /i, • • • , /at-i- 
Since 

as |A| — > co, we have for each n e {0, 1, . . . , N — 1}, 



i-/ r {E«„ m (A,/ m },A + / r O(J 5 ),A, 



/» = / 

where r is a sufficiently large positive number, T r is the circle in the A-plane of radius r centered 
at the origin. 

Denote by Ai, . . . , X p all the distinct roots of the polynomial Pjv(A) (which coincides by (|2.6p 
with the characteristic polynomial of the matrix J up to a constant factor) and by mi, . . . , m p 
their multiplicities, respectively: 

P N (X) = c(A - Ai) mi • • • (A - A p )"\ (5.8) 
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where c is a constant. We have 1 < p < N and mi + • • • + m p = N. By (j5.8|) . we can rewrite the 
rational function (3at(A)/-P/v(A) as the sum of partial fractions: 

where /3kj are some uniquely determined complex numbers depending on the matrix J. Substi- 
tuting (15. 6h in (|5.7[) and taking into account (|5.5I) . (|5.9I) we get, applying the residue theorem 
and passing then to the limit as r — ► oo, 

^ = EEr%T^FT[ F ( A ) P «( A )] > n€{0,l,...,iV-l}, (5.10) 

fe=l j=l V J- L J A=A fc 

where 

AT-l 

^(A) = 2 /mi'm(A). (5.11) 
m=0 

Now define on C2iv[A] the functional f2 by the formula 

fe=ij=i U >' 
Then formula (|5.10p can be written in the form 

f n = (n,F(X)P n (X)), n€{0,l,...,JV-l}. (5.13) 
From here by (|5,lip and the arbitrariness of {/ m }^ZQ it follows that the "orthogonality" relation 

(n,P m (X)P n (X)) = 5 mn , m,ne{0,l,...,JV-l}, (5.14) 
holds. Further, in virtue of (|5.8p and f|5. 12j) we have also 

(n,p m (\)p N (\)) =0, me{0,i,...,jv}. (5.15) 

These mean by Theorem [1] that the generalized spectral function of the matrix J has the 
form (I5TT2D . ■ 

Definition 2. The collection of the quantities 

{Afc,/3fcj (j = 1,... ,m k ,k = 1, . . . 

determining the structure of the generalized spectral function of the matrix J according to 
Theorem HI we call the spectral data of the matrix J. For each k G {1, . . . ,p} the sequence 

{Pkli • • • i Pkm k } 

we call the normalizing chain (of the matrix J) associated with the eigenvalue (the sense of 
"normalizing" will be clear below in Section [8]). 
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If we delete the first row and the first column of the matrix J given in (jl.ip , then we get the 
new matrix 



6 (1) a {1) 







a« b? a« 



a 



(i) ii) 

1 °2 















. 


• • °N-4 


a (1) 











• 


•• a (1) 




a {1) 
TV— 








• 





a (1) 

N—3 


b {1) 
N— 



where 



a (1) - a ,1 



n£ {0,l,...,iV-3}, 
n€ {0,l,...,JV-2}. 



The matrix jW is called the /irsi truncated matrix (with respect to the matrix J). 

Theorem 5. The normalizing numbers fikj of the matrix J can be calculated by decomposing 
the rational function 

det(jW - XI) 
det( J - XI) 

into partial fractions. 

Proof. Let us denote the polynomials of the first and the second kinds, corresponding to the 
matrix by P^\X) and Qn\X), respectively. It is easily seen that 



P« (A) = aoQn+iW, n € {0, 1, . . . , N - 1}, 

Q { n ] W = -{(A - b )Q n+1 (X) - P n+1 (X)}, 
ao 



ne {0,1,... ,N- 1}. 

Indeed, both sides of each of these equalities are solutions of the same difference equation 
<4-i?/n-i + b$y n + a^Vn+i = Xy n , n G {0, 1, . . . , N - 2}, 



(5.16) 
(5.17) 



a (1) -1 



and the sides coincide for n = — 1 and n = 0. Therefore the equality holds by the uniqueness 
theorem for solutions. 

Consequently, taking into account Lemma Q] and using (|5.16p , we have 

det (j( 1 ) - XI) = (-If^a^a? ■ ■ • aJJ } _ 3 ^v-i(A) = (-1)^1 ■ ■ ■ a N _ 2 a Q N (X). 

Comparing this with (|2.6p . we get 

Qn(X) _ det(j( 1 ) - XI) 
P N (X) ~ det( J - XI) 

so that the statement of the theorem follows from (15.91). ■ 
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6 Inverse problem from the spectral data 

By the inverse spectral problem is meant the problem of recovering matrix J, i.e. its entries a n 
and b n , from the spectral data. 

Theorem 6. Let an arbitrary collection of complex numbers 

{Afc,/3fcj (i = 1, ... ,nik,k = 1, ... ,p)} (6.1) 

be given, where Ai, A2, . . . , X p (1 < p < N) are distinct, 1 < < N , and m\ + • • • + m p = N . 
In order for this collection to be the spectral data for some Jacobi matrix J of the form (|1.1|) 
with entries belonging to the class (|1.2p . it is necessary and sufficient that the following two 
conditions be satisfied: 

(*) £&i = i; 

k=l 

(ii) D n 7^ 0, for n € {1, 2, . . . , N — 1}, and Dn = 0, where D n is defined by (|3.24p in which 
EE (6-2) 



si = 

k=i j=i 



min{mfc,Z + 1} ; ( is a binomial coefficient. 



Proof. The necessity of conditions of the theorem follows from Theorem [3] because the generali- 
zed spectral function of the matrix J is defined by the spectral data according to formula (15. ip 
and therefore the quantity (|6.2p coincides with ^fi, A^) . Besides, 

p 

Y^Pki = (n,l) = s = D . 

k=l 

Note that the condition (ra) of Theorem [2] holds with 

T(A) = (A- Ai) mi •••(A- X p ) mp . (6.3) 



Let us prove the sufficiency. Assume that we have a collection of quantities (|6.ip satisfy- 
ing the conditions of the theorem. Using these data we construct the functional 0, on C2jv[A] 
by formula (|5.ip . Then this functional satisfies the conditions of Theorem [3] and therefore 
there exists a matrix J of the form (II. ip . (11. 2p for which is the generalized spectral function. 
Now we have to prove that the collection (|6.ip is the spectral data for the recovered matrix J. 
For this purpose we define the polynomials P_i(A), -Po(A), • • • ,P/v(A) as the solution of equa- 
tion (|5,2p . constructed by means of the matrix J, under the initial conditions (|5.3p . Then the 
relations (|2.7p . (|2.8p and the equalities 



o n = (fi,AP n (A)P n+ i(A)), n€{0,l,-..,JV-2}, (6.4) 
b n = (n,\P^(X)), n€{0,l,-..,JV-l} (6.5) 



hold. We show that (|5.8p holds which will mean, in particular, that Ai, . . . , A p are eigenvalues 
of the matrix J with the multiplicities mi, . . . , m p , respectively. 

Let T(A) be defined by (|6.3[) . Let us show that there exists a constant c such that 

av- 2 Pjv- 2 (A) + &tv-iPtv-i(A) + cT(A) = XP N -i(X) (6.6) 

for all A G C. If we prove this, then from here and (|5.2p with y^ = Pfc(A) and n = N — 1 we get 
that Pjv(A) = cT(A). 
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Since degP n (A) = ra(0<n<iV — 1), degT(A) = mi + • • • + m p = N, the polynomials 
Po(A), . . . , Pjv-i(A), T(A) form a basis of the linear space of all polynomials of degree < N. 
Therefore we have the decomposition 

N-l 

APv-i(A) = cT(A) + ^ Cn P n (A), (6.7) 

n=0 

where c, co,c\, . . . , cn-i are some constants. By (|6.3p and (|5.1|) it follows that 

(Q,T(X)P n (X)) = 0, n€{0,l,...,JV}- 
Hence taking into account the relations (|2.7p . ()2.8[) and (|6.4p . (|6.5p . we find from (|6.7p that 

C n = (0<?l<iV-3), CAr_ 2 =aAr_2, Cjv_i = 6jv-1- 

So (|6.6p is shown. 

It remains to show that for each k G {1, . . . ,p} the sequence {/?&!, . . . , f3km k } is the normali- 
zing chain of the matrix J associated with the eigenvalue A&. Since we have already shown that 
Afc is an eigenvalue of the matrix J of the multiplicity m^, the normalizing chain of J associated 
with the eigenvalue A& has the form {(3ki, ■ ■ ■ , /3fcm fc }- Therefore for (Q, G(A)) we have an equality 
of the form (|5,ip in which fl^j is replaced by (3^. Subtracting these two equalities for (Q,G(X)} 
each from other we get that 

tETrf G(H) ^) = ° fOTa11 G(X)eC 2N [X}. 
k=ij=i ^ >' 

bince the values G^-^(X k ) can be arbitrary numbers, we get that fikj = Pkj for all k and j. ■ 

Under the conditions of Theorem [6] the entries a n and b n of the matrix J for which the 
collection (I6.ip is spectral data, are recovered by formulas (I3.35p . ()3.36p . 

7 Characterization of generalized spectral functions 
of real Jacobi matrices 

In this section, we characterize generalized spectral functions of real Jacobi matrices among 
the generalized spectral functions of complex Jacobi matrices. Let m be a nonnegative integer. 
Denote by M2m[A] the ring of all polynomials in A of degree < 2m with real coefficients. 

Definition 3. A linear functional defined on the space C2m[A] is said to be positive if 

(n,G(X)) > o 

for all polynomials G(X) £ M2m[A], which are not identically zero and which satisfy the inequality 

G(X) > 0, -oo < A < oo. 

Lemma 2. If Q is a positive functional on C2 m [A], then it takes only real values on lR2m[^]- 

Proof. Since the functional f2 is positive, the values (fi, A 2fc ), k G {0, 1, . . . , m} are real (more- 
over they are positive). Next, the monomial A 2fc_1 , k G {1, 2, . . . , m} is represented as a difference 
of two nonnegative polynomials of degree 2k : 

Therefore the values ($1, A 2fc_1 ), k G {1, 2, . . . , m} are also real to be a difference of two positive 
numbers. Thus, (fi, A n ) is real for any n G {0, 1, . . . , 2m}. Hence (f2,G(A)} is real for any 
G(X) G M 2m [A]. ■ 
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Lemma 3. A linear functional O on C2 m [A] is positive if and only if D n > for all n G 
{0, 1, . . . , m}, where 



in which 



s Si 
si s 2 

Sn 



Sn 
Sn+1 

S2n 



n = 0, 1, . . . , m, 



si = (n,X l ), l = 0,1,..., 2m. 

Proof. Any polynomial G(X) 6 M2m[A] which is not identically zero and which satisfies the 
inequality 



G(X) > 0, -oo < A < oo, 
can be represented in the form 
G(A) = [A(X)] 2 + [5(A)] 2 , 



(7.1) 



(7.2) 



where ^4(A), B(X) are polynomials of degrees < m with real coefficients. Indeed, it follows 
from ()7.ip that the polynomial G(A) has even degree: degG(A) = 2p, where p < m. Therefore 
its decomposition into linear factors has the form 

p 

G(X) = cJJ(A - a k - i(3 k )(X - a k + i/3 k ), 

k=l 

where c > 0, (3 k > 0, a k are real (among the roots a k + ifiki °f course, may be equal). Now 
setting 

v 

^Tc\{{X-a k -ih)=A{X)+iB{X), 

k=l 

we get that the polynomials A(X), B(X) have real coefficients and ()T.2f) holds. 
Now writing 

p p 
A(X) = ^2x k X k , B(X) = y £y k X k , 

k=l k=l 

where x k , y k are real numbers, we find 

p p 
(Q,G(X)} = ^2 Sj+kXjX k + ^2 Sj+kVjVk- 

j,k=0 j,k=0 

This implies the statement of the lemma. ■ 

Theorem 7. In order for a given linear functional Q, on C-2n[X], to be the generalized spectral 
function for a real Jacobi matrix of the form (jl.ip . (|1.3|) it is necessary and sufficient that the 
following three conditions be satisfied: 

(i) <fi, 1) = 1; 

(ii) Q is positive on C2JV-2M; 
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(Hi) there exists a polynomial T(X) of degree N such that (Q, G(X)T(X)) = for all polynomials 
G(X) with degG(A) < N. 

Proof. Necessity. The condition (£1, 1) = 1 follows from (|2.7|) with m = n = 0. To prove 
positivity on C2at-2[A] of the generalized spectral function Q of the real Jacobi matrix J, take 
an arbitrary polynomial G(X) E 1^2^-2 [A] which is not identically zero and which satisfies the 
inequality 

G(X) > 0, — oo < A < co. 

This polynomial can be represented in the form (see the proof of Lemma [3]) 

G(X) = [A(X)] 2 + [5(A)] 2 , (7.3) 

where ^4(A), B(X) are polynomials of degrees < N — 1 with real coefficients. Since the poly- 
nomials -Po(A), A (A), . . . ,P/v_i(A) have real coefficients (because J is a real matrix) and they 
form a basis of Mtv_i[A], we can write the decompositions 

N-l N-l 

A(X) = ]T c k P k (X), B(X) = d k P k (X), 

k=l k=l 

where c k , d k are real numbers not all zero. Therefore using the "orthogonality" property (12. 7ft 
we get from (|7.3p . 



N-l 

(n,G(X))= J2(4 + 4)>o. 

j,k=0 

The property of f2 indicated in the condition (Hi) of the theorem follows from (|2.8p if we take 
T(X) = P N (X). 

Sufficiency. It follows from the conditions of the theorem that all the conditions of Theorem [2] 
are satisfied. In fact, we need to verify only the condition (ii) of Theorem (2) Let for some 
polynomial G(A), degG(A) = n < N — 1, 

(fi, G(X)H(X)) = (7.4) 

for all polynomials H(X), degH(X) = n. We have to show that then G(X) = 0. Setting 

n n 

G(X) = Y,9kX k , H(X) = Y,h j V, 

k=0 j=0 

we get from (|7.4p that 

n / n \ 

j=0 \k=0 / 

Since ho, h±, . . . , h n (h n ^ 0) are arbitrary, the last equation gives 

n 

^9kSk+j = 0, j = 0, l,...,n. (7.5) 

k=0 

This is a linear homogeneous system of algebraic equations with respect to go, g±, . . . , g n and the 
determinant of this system coincides with the determinant D n . Prom the condition (ii) of the 
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theorem it follows by Lemma [3] that D n > 0. So D n ^ and hence system (j7.5|) has only the 
trivial solution go = gi = • • • = g n - 

Thus, all the conditions of Theorem[2]are satisfied. Therefore there exists, generally speaking, 
a complex Jacobi matrix J of the form (jl.ip . (jl.2p for which ft is the generalized spectral 
function. This matrix J is constructed by using formulas (|3,35p . (|3.36[) . It remains to show that 
the matrix J is real. But this follows from the fact that by Lemma [2] and Lemma [3] we have 
D n > for n £ {0, 1, . . . , N — 1} and the determinants A n are real. Therefore formulas f|3.35j) . 
(13. 36ft imply that the matrix J is real. ■ 

If we take into account Lemma [3j then it is easily seen from the proof of Theorem [3] that 
Theorem [7] is equivalent to the following theorem. 

Theorem 8. In order for a given linear functional f2, defined on C2n[X], to be the generalized 
spectral function for some real Jacobi matrix J of the form with entries belonging to the 

class (|1.3p . it is necessary and sufficient that 

D = l, D n >0 (n = l,2,...,JV-l), and D N = 0, 

where D n is defined by (|3.24p and (|3.23p . 

Under the conditions of Theorem [8] the entries a n and b n of the matrix J for which the 
functional is the spectral function are recovered by formulas (|3.35p . (|3,36p , 

8 Structure of generalized spectral functions 
of real Jacobi matrices 

First we prove two lemmas which hold for any complex Jacobi matrix J of the form (jl.ip , (|1.2p . 
Having the matrix J consider the difference equation (15. 2p and let {P n (X)}^ = _ 1 and {Q n (X)}^ = _ 1 
be the solutions of this equation satisfying the initial conditions (|5.3p and (|5.4p . respectively. 



Lemma 4. The equation 

P JV _i(A)Q 7V (A) - PnWQn-iW = 1 (8.1) 

holds. 

Proof. Multiply first of the equations 

a n _iP n _i(A) + b n P n {\) + a n P n+l {\) = AP n (A), (8.2) 
n G {0, 1, . . . , N — 1}, a_i = a^-i = 1, 
On-lQn-l(A) + b n Q n (X) + a n <5n+i(-^) = AQ n (A), 
77, 6 {0, 1, . . . , N — 1}, a_i = oat_i = 1, 

by Qn(A) and second by P n (A) and subtract the second result from the first to get 

a n _![P n _i(A)Q n (A) - P„(A)Q„_i(A)] 

= a n [P n (A)Q n+1 (A) - P n+ i(A)Q n (A)], n G {0, 1, . . . , N - 1}. 

This means that the expression (Wronskian of the solutions P n (A) and Q n (\)) 

a n [P n (X)Q n+1 (X) - P n+1 (X)Q n (X)] 

does not depend on n € {—1,0, 1, . . . , N — 1}. On the other hand the value of this expression at 
n = — 1 is equal to 1 by (|5,3p . (|5.4p . and a_i = 1. Therefore 

a n [P n (X)Q n+ i(X) - P n+1 (X)Q n (X)} = 1 for all n G {-1, 0, 1, . . . , N - 1}. 

Putting here, in particular, n = N — 1, we arrive at (|8,ip . ■ 
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Lemma 5. The equation 

JV-i 



P N ^(\)P' N (X) - P N (X)P' N _ 1 (X) = P nW (8-3) 

n=0 

holds, where the prime denotes the derivative with respect to X. 
Proof. Differentiating equation (|8.2p with respect to A, we get 

a n - 1 P' n _ 1 (X) + KKW + «nP;+i(A) = AP^(A) + P n (A), (8.4) 
n 6 {0, 1, . . . ,N - 1}, a_i = oat-i = 1. 

Multiplying equation (|8.2|) by P^(A) and equation (|8.4|) by P n (A), and subtracting the left and 
right members of the resulting equations, we get 

a n _ 1 [P n _ 1 (A)i>UA) - P n (A)P,U(A)] - a n [P n (X)P' n+1 (X) - P n+1 (A)P^(A)] = -P n 2 (A), 
n £ {0,1,..., N- 1}. 

Summing the last equation for the values n = 0,1, ... ,m (m < N — 1) and using the initial 
conditions (15.31). we obtain 



a m [P m (A)P^ +1 (A) - P m+1 (X)P^(X)} = ^n(A), me{0,l,...,N-l}. 

n=0 

Setting here, in particular, m = N — 1 and taking into account ajy-i = 1, we get (|8.3p . ■ 
Now we consider real Jacobi matrices of the form (jl.ip . (jl.3p . 

Lemma 6. For any real Jacobi matrix J of the form (jl.ip . (jl.3p i/ie roots o/ i/ie polyno- 
mial P^ (A) are simple. 

Proof. Let Ao be a root of the polynomial P/v(A). The root Ao is an eigenvalue of the matrix J 
by ([2.6p and hence it is real by Hermiticity of J. Putting A = Ao in (|8.3p and using P/v (Ao) = 0, 
we get 

AT-l 

P JV _ 1 (Ao)P^(A ) = Y P n 2 (A ). (8.5) 

n=0 

The right-hand side of (|8.5p is different from zero because the polynomials P n (A) have real 
coefficients and hence are real for real values of A, and Po(A) = 1. Consequently P^(Ao) 7^ 0, 
that is, the root Ao of the polynomial P/v(A) is simple. ■ 

Lemma 7. Any real Jacobi matrix J of the form (jl.lj) . (I1.3P /ms precisely N real and distinct 
eigenvalues. 

Proof. The reality of eigenvalues of J follows from its Hermiticity. Next, the eigenvalues of J 
coincide, by (|2.6p . with the roots of the polynomial Pat (A). This polynomial of degree N has N 
roots. These roots are pairwise distinct by Lemma [6j ■ 

The following theorem describes the structure of generalized spectral functions of real Jacobi 
matrices. 
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Theorem 9. Let J be a real Jacobi matrix of the form (jl.ip . (| 1 .3 j) and Q be its generalized 
spectral function. Then for any G(X) 6 C2jv[A] 



iV 



(n,G(X)) = Y,PkG(\ k ), (8.6) 
fe=i 

where X\, . . ., Aj\r are £/ie eigenvalues of the matrix J and /3\, . . . ,{3n <ire positive real numbers 
uniquely determined by the matrix J. 

Proof. By Lemma [6l the roots Ai, . . . , Ajy of the polynomial P/v(A) are simple. Therefore the 
formula (|8.6j) follows from (|5.1h and the decomposition (|5.9p takes the form 

Qn(X) _ (3k 
p nW A-Afc' 

Hence 

Qjv(A fc ) =/3 fc P^(A fc ). (8.7) 

On the other hand, putting A = A& in (|8.ip and (|8.3p and taking into account that -Pa^A/J = 0, 
we get 



J P7v-i(A fe )Q7v(A fc ) = 1, 



N-l 



P N -i(X k )P^(X k ) = £ P„ 2 (A fc ), (8.9) 
n=0 

respectively. Comparing (|8.7p . (|8.8p . and (|8.9p . we find that 

^ = {E P -( A ^} . ( 8 - 10 ) 

ln=0 J 

whence we get, in particular, that (3k > 0. ■ 

Since {P n (A/ c )}^ r =0 1 is an eigenvector of the matrix J corresponding to the eigenvalue X k , it 
is natural, according to the formula (|8.10p . to call the numbers (3 k the normalizing numbers of 
the matrix J. 

Definition 4. The collection of the eigenvalues and normalizing numbers 

{X k ,(3 k (k = l,...,N)} 
of the matrix J of the form (jl.ip . (jl.3|) we call the spectral data of this matrix. 

Remark 2. Assuming that Ai < A2 < • • • < A v, let us introduce the nondecreasing step 
function uj(X) on (—00, 00) by 

LO(X) = ^, 
\ k <\ 

where u>(X) = if there is no A^ < A. So the eigenvalues of the matrix J are the points of 
increase of the function uj(X). Then equality (18. 6|) can be written as 



(0,G(A)) = / G(X)du)(X), 
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where the integral is a Stieltjes integral. Therefore the orthogonality relation (|5,14p can be 
written as 

/oo 
P m {X)P n {X)dw{X) = 5 mn , m,7ie{0,l,...,j\T-l} 
-oo 

and the expansion formula f)5 . 13|) as 

/oo 
F(X)P n (X)dco(X), n€{0,l l ...,JV-l} J 
-oo 

where F(X) is defined by (|5.11|) : 

N-l 

Such function w(A) is known as a spectral function (see, e.g., [21]) of the operator (matrix) J. 
This explains the source of the term "generalized spectral function" used in the complex case. 



9 Inverse spectral problem for real Jacobi matrices 

By the inverse spectral problem for real Jacobi matrices we mean the problem of recovering the 
matrix, i.e. its entries, from the spectral data. 

Theorem 10. Let an arbitrary collection of numbers 

{X k ,f3 k (k = l,...,N)} (9.1) 

be given. In order for this collection to be the spectral data for a real Jacobi matrix J of the 
form (|l.ip with entries belonging to the class (|1.3p . it is necessary and sufficient that the following 
two conditions be satisfied: 

(i) The numbers X±, . . . , A/v are real and distinct. 

N 

(ii) The numbers @\,... ,0n are positive and such that Pk = 1- 

k=l 

Proof. The necessity of the conditions of the theorem was proved above. To prove the sufficien- 
cy, assume that we have a collection of quantities (|9.ip satisfying the conditions of the theorem. 
Using these data we construct the functional 0, on C2tv[A] by the formula 

TV 

(n, G(A)> = £ (3 k G(X k ), G(X) e C 2N [X}. (9.2) 

k=l 

Then this functional O satisfies the conditions of Theorem [71 Indeed, we have 

/v 



k=l 



Next, let G(A) G M2/v-2[A] be an arbitrary polynomial which is not identically zero and which 
satisfies the inequality 

G(X) > 0, -oo < A < oo. 
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This polynomial can be represented in the form (see the proof of Lemma [3]) 

G(A) = [A(X)f + [5(A)] 2 , 
where -A(A), B(X) are polynomials of degrees < N — 1 with real coefficients. Then 

TV N TV 

<fi,G(A)> = ^/3 fc G(A fe ) = ^/3 fe [^(A fe )] 2 + ^/3,[£(A fe )] 2 > 0. (9.3) 

k=l k=l k=l 

We have to show that the equality sign in (|9.3p is impossible. If we have the equality in (|9.3p . 
then, since all the are positive, we get 

A(Ai) = • • • = A(Ajv) = and B(Ai) = • • • = S(Atv) = 0. 

Hence A(X) = and B(X) = because Ai, . . . , Aat are distinct and deg^4(A) < N—l, degB(X) < 
N — 1. Therefore we get G(X) = which is a contradiction. Finally, if we take 

T(A) = (A-A 1 )---(A-A 7V ), 
then the condition (Hi) of Theorem [7] is also satisfied: for any polynomial G(X), 

TV 

(fi, G(A)T(A)> = ^ /3 fc G(A fc )T(A fc ) = 0. 
fc=i 

Thus, the functional f2 defined by the formula (19. 2D satisfies all the conditions of Theorem [71 
Therefore there exists a real Jacobi matrix J of the form (jl.ip . (| 1 .3 j) for which f2 is the generalized 
spectral function. Further, from the proof of sufficiency of the conditions of Theorem[6]it follows 
that the collection {A^, (3^ {k = 1, . . . , N)} is the spectral data for the recovered matrix J. ■ 

Note that under the conditions of Theorem [10] the entries a n and b n of the matrix J for which 
the collection ()9. 1 j) is spectral data, are recovered by formulas (|3.35p . (|3.36p . 

Acknowledgements 

This work was supported by Grant 106T549 from the Scientific and Technological Research 
Council of Turkey (TUBITAK). 

References 

[1] Boley D., Golub G.H., A survey of matrix inverse eigenvalue problems, Inverse Problems 3 (1987), 595-622. 

[2] Ikramov Kh.D., Chugunov V.N., Inverse matrix eigenvalue problems, J. Math. Sciences 98 (2000), 51-136. 

[3] Chu M.T., Golub G.H., Inverse eigenvalue problems: theory, algorithms, and applications, Oxford University 
Press, New York, 2005. 

[4] Marchenko V.A., Expansion in eigenfunctions of non- self adjoint singular second order differential operators, 
Mat. Sb. 52 (1960), 739-788 (in Russian). 

[5] Rofe-Beketov F.S., Expansion in eigenfunctions of infinite systems of differential equations in the non- 
selfadjoint and selfadjoint cases, Mat. Sb. 51 (1960), 293-342 (in Russian). 

[6] Guseinov G.Sh., Determination of an infinite non-selfadjoint Jacobi matrix from its generalized spectral 
function, Mat. Zametki 23 (1978), 237-248 (English transl.: Math. Notes 23 (1978), 130-136). 

[7] Guseinov G.Sh., The inverse problem from the generalized spectral matrix for a second order non-selfadjoint 
difference equation on the axis, Izv. Akad. Nauk Azerb. SSR Ser. Fiz.-Tekhn. Mat. Nauk (1978), no. 5, 16-22 
(in Russian). 



28 



G.Sh. Guseinov 



Kishakevich Yu.L., Spectral function of Marchenko type for a difference operator of an even order, Mat. 
Zametki 11 (1972), 437-446 (English transl.: Math. Notes 11 (1972), 266-271). 

Kishakevich Yu.L., On an inverse problem for non-selfadjoint difference operators, Mat. Zametki 11 (1972), 
661-668 (English transl.: Math. Notes 11 (1972), 402-406). 

Bender CM., Making sense of non-Hermitian Hamiltonians, Rep. Progr. Phys. 70 (2007), 947 1018, 
|hep-th/0703096l 

Znojil M., Matching method and exact solvability of discrete PT-symmetric square wells, J. Phys. A: Math. 
Gen. 39 (2006), 10247-10261, |quant-ph/0605209| 

Znojil M., Maximal couplings in PT-symmetric chain models with the real spectrum of energies, J. Phys. A: 
Math. Theor. 40 (2007), 4863-4875, |math-ph/0703070| 



Znojil M., Tridiagonal PT-symmetric N by N Hamiltonians and fine-tuning of their observability domains 
in the strongly non-Hermitian regime, J. Phys. A: Math. Theor. 40 (2007), 13131-13148. larXiv:0709. 15691 

Allakhverdiev B.P., Guseinov G.Sh., On the spectral theory of dissipative difference operators of second 
order, Mat. Sb. 180 (1989), 101-118 (English transl.: Math. USSR Sbornik 66 (1990), 107-125). 

Guseinov G.Sh., Completeness of the eigenvectors of a dissipative second order difference operator, J. Dif- 
ference Equ. Appl. 8 (2002), 321-331. 

van Moerbeke P., Mumford D., The spectrum of difference operators and algebraic curves, Acta Math. 143 
(1979), 93-154. 

Sansuc J.J., Tkachenko V., Spectral parametrization of non-selfadjoint Hill's operators, J. Differential Equa- 
tions 125 (1996), 366-384. 

Egorova I., Golinskii L., Discrete spectrum for complex perturbations of periodic Jacobi matrices, J. Dif- 
ference Equ. Appl. 11 (2005), 1185-1203, rmath.SP/0503627| 

Atkinson F.V., Discrete and continuous boundary problems, Academic Press, New York, 1964. 

Akhiezer N.I., The classical moment problem and some related questions in analysis, Hafner, New York, 
1965. 

Berezanskii Yu.M., Expansion in eigenfunctions of selfad joint operators, Translations of Mathematical Mono- 
qraphs, Vol. 17, American Mathematical Society, Providence, R.I., 1968. 

Nikishin E.M., Sorokin V.N., Rational approximations and orthogonality, Translations of Mathematical 
Monoqraphs, Vol. 92, American Mathematical Society, Providence, R.I., 1991. 

Teschl G., Jacobi operators and completely integrable nonlinear lattices, Mathematical Surveys and Mono- 
qraphs, Vol. 72, American Mathematical Society, Providence, R.I., 2000. 



